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ABSTRACT 



We study the response of the low-degree, solar p-mode frequencies to the unusually extended minimum of solar surface activity since 
2007. A total of 4768 days of observations collected by the space-based, Sun-as-a-star helioseismic GOLF instrument are analyzed. 
A multi-step iterative maximum-likelihood fitting method is applied to subseries of 365 days and 91.25 days to extract the p-mode 
parameters. Temporal variations of the / = 0, 1, and 2 p-mode frequencies are then obtained from April 1996 to May 2009. While the 
p-mode frequency shifts are closely correlated with solar surface activity proxies during the past solar cycles, the frequency shifts of 
the / = and / = 2 modes show an increase from the second half of 2007, when no significant surface activity is observable. On the 
other hand, the 1=1 modes follow the general decreasing trend of the solar surface activity. The different behaviours between the 
1 = and 1 = 2 modes and the / = 1 modes can be interpreted as different geometrical responses to the spatial distribution of the solar 
magnetic field beneath the surface of the Sun. The analysis of the low-degree, solar p-mode frequency shifts indicates that the solar 
activity cycle 24 started late 2007, despite the absence of activity on the solar surface. 
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1. Introduction 

The success of helioseismology is due to its capability to accu- 
rately measure the p-mode parameters of the solar eigenmode 
spectrum, which allow us to infer unique information about the 
internal structure and dynamics of the Sun from its surface all 
the way down to the core. It has greatly contributed to a better 
understanding of the Sun and provided insights into the com- 
plex solar magnetism, through for instance the variability of the 
characteristics of the p-mode spectrum. 

Evidence of low-degree p-mode frequency chang es with so- 
lar activity, first revealed bvlWoodard & Novesl d 1985b . were es- 
tablished bv lPalle et al.l d 1989b with the analysis of helioseismic 
observati ons spanning the complete sol ar cycle 21 (1977-1988). 
Later on, Woo dard & Libbrechtl dl99ll) measured similar varia- 
tions in intermediate and high angular-degree modes. As longer, 
higher quality, and continuous helioseismic observations became 
available, the solar p-mode frequencies proved to be very sensi- 
tive to the solar surface activity not only during solar cycle 23 
(see for example. iHowe et^l2002t|jim enez-Re ves et al.Ll2004 
and references therein) but also over the last thr ee solar cycles 
(21, 2 2, and 23), spanning a period of thirty years dChaplin et al.L 
120071) . with high levels of correlation with solar surface activity 
proxies. 

The frequency shifts have also been show n to be angular- 
degree dependent dJimenez-Reves et"ail 1200 lb . or rather mode- 
inertia dependent, but IHowe et al.l (2002) found that the latitu- 
dinal distribution of the frequency differences for high / modes 
shows actually close temporal and spatial correlation with the 
surface magnetic field distribution. A marginal but significant / 
dependence of the frequency shifts of the low-degree p modes 



has also been clear ly uncovered dJimenez-Reves et"aTl [2004; 
ChaphnetaU 12001 . The size of the shift also reflects differ- 
ences in mode inertia. However, the change in inertia with I is 
small for low-degree modes, and the spatial distribution domi- 
nates as well. 

The unusual long and deep solar activity minimum as well 
as the lack of a visible onset of the new cycle 24 have motivated 
the present work in order to investigate the response of the low- 
degree p-mode frequencies to this extended minimum (2007 till 
present) as compared with the rest of cycle 23 (since 1996). 



2. Data and analysis 

2.1. Data sets 

Observations collected by the space-based instrument Global 
Oscillations at Low Frequency (GOLF) instrument onboard the 
Solar and Heliospheric Observatory (SOHO) spacecraft were 
used for this work. GOLF is a resonant scattering spectropho- 
tometer measuring the Doppler wavelength shift - integrated 
over the solar surface - in the Di and D 2 Fraunhofer sodium 
lines at 589.6 and 589.0 nm respectively dGabriel et all Il995l 
A total of 4768-day velocity time series dGarcfa et all [2005 
st arting on 1996 April 11 w ere used (see calibration method 
in Jimenez ~Reyes et a"Ql2003l) . Contiguous subseries of 365 and 
91.25 days were produced with a four-time overlap. The series 
with a filling factor less than 60% for the 365-day series and 70% 
for the 91.25-day series were not used in the following analysis, 
resulting in a total of 49 and 197 non-independent time series of 
365 days and 91.25 days respectively. Thus, the averaged duty 
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cycles for the 365-day and 91.25-day time series were 94.9% 
and 97.8% respectively. 

We have also considered contemporaneous observations 
collected by the grou nd-based Mark-I in strument at the 
Observatorio del Teide dBrookes et all Il978l) . Mark-I is also 
a resonant scattering spectrophotometer but observing in the 
potassium Fraunhofer line at 769.9 nm. A total of 19 indepen- 
dent yearly time series of the best consecutive 108 days (mean 
duty cycle 38%) were used. 



2.2. Mode parameter extraction 

The power spectrum of each 365-day and 91.25-day time se- 
ries was fitted to yield estimates of the mode parameters. 
This fitting was performed by a mul ti-step iterative method 
dGellv et all 120021; ISalabert et al.L 120071) . The asymmetric pro- 
file of Nigam & Kosovichevl ( 19981) was used to describe each 
component, as: 



P nJ (x) = H n , 



(1 + a„jx nJ ) 2 + a 2 nl 
1 + x 2 , 

n,l 



+ B, 



(1) 



where x n j = 2(v - v„j)/r„j, and v n j, F n j, and H n j represent 
the central frequency, the Full-Width-at-Half-Maximum (fwhm), 
and the power height of the spectral density respectively for a 
given mode (n, /). The peak asymmetry is described by the pa- 
rameter a n ,i, while B represents an additive, constant background 
level in the fitted window. Due to their close proximity in fre- 
quency, modes are fitted in pairs (i.e., / = 2 with 0, and I = 3 
with 1). While each mode parameter within a pair of modes are 
free, the peak asymmetry is set to be the same within pairs of 
modes. The mode parameters were extracted by maximizing a 
likelihood function following the x 2 with 2 degrees-of-freedom 
statistics of the power spectrum. The natural logarithms of the 
mode height, fwhm, and background noise were varied which 
result in normal distributions. The formal uncertainties on each 
parameter were then derived from the inverse Hessian matrix. A 
similar analysis was applied to the Mark-I spectra. 

2.3. Determination of the frequency shifts 

The temporal variations of the p-mode frequencies were ob- 
tained by comparing each measured (I, n) frequency with the cor- 
responding (/, n) averaged frequency over the whole set of ana- 
lyzed spectra for the different datasets. The formal uncertainties 
returned by the fits were used as weights in the average computa- 
tion. The frequency shifts thus obtained were then averaged over 
the frequency range 2000 ^Hz < v < 3300 /iHz. Note that due to 
its lower signal-to-noise ratio, we did not use the I = 3 modes in 
the following analysis. Linear regressions between the frequency 
shifts and solar surface activity proxies were performed. Mean 
values of daily measurements of the 10.7-cm radio flux and the 
Mount Wilson Magnetic Plage Strength Index (MPSI) were cal- 
culated over the same 365-day and 91.25-day subseries. 



3. Results 

As the minimum of solar cycle 23 has been longer than predicted 
and that no surface activity is being observed since the end of 
2007, we might expect that the frequency shifts would also show 
an extended minimum. Figure Q] shows the temporal variation of 
the frequency shifts (<5vo i 2) averaged over the modes / = 0, 1, 



and 2 extracted from the analysis of both the 365-day and 91.25- 
day GOLF spectra. As already widely reported, the frequency 
shifts represented on Fig. Q] are well correlated with the solar 
surface activity over cycle 23 (1996-2007). However, as we can 
see on Fig. [1] the mean frequency shifts start increasing from 
the second half of 2007 while no significant surface activity was 
observed. Excellent qualitative agreement was obtained with the 
Mark-I data, detailed work being in progress. 

Figure [2] shows the dependence of the frequency shifts (6vi) 
on the angular degree I with time, measured from the 365-day 
spectra. The frequency shifts at each / per unit of change in two 
global solar activity proxies, as well as the associated Pearson 
correlation coefficient, r p , the Spearman rank correlation coef- 
ficient, r s , and the two-sided significance of its deviation from 
zero, P s are given in Table Q] These values were obtained from 
the analysis of the 91.25-day spectra by using only independent 
points (43 for the period April 1996-June 2007, and 7 for the 
period June 2007-May 2009). Note that they correspond to one 
combination of independent series and that they are consistent 
within lcr with any other chosen sets. Given the range of ana- 
lyzed frequencies and the considered time spans, the amount of 
the uncovered shifts over cycle 23 a nd their correlations are con- 
sist ent within the uncerta inties with Jimene z~Reves et ail ([2004) 
and lChaplin et all (|2004). The / = 2 modes have larger variations 
in comparison to the I = 1 and I = modes respectively. 

From the end of 2007, the frequency shifts of the I = 1 
modes keep slowly decreasing, following the trend of the so- 
lar surface activity, while the I = and I = 2 modes have a 
very different behaviour. Indeed, the modes I = and 2 show 
increasing frequency shifts, such as expected if solar cycle 24 
has already started to progress. The / = and 2 frequency shifts 
show negative correlation with surface activity, while the I = 1 
frequency shifts are in some extent still correlated with surface 
activity (Table[T|i. Note that the significances of having null cor- 
relation are larger because of much fewer number of independent 
points (7 points) for the period June 2007-May 2009. In order 
to stress this different behaviour with /, Fig.[3]compares the fre- 
quency shifts of the / = and I = 1 modes in relation with those 
of the I - 2 modes for the two periods April 1996-June 2007 
and June 2007-May 2009, confirming what is shown on Fig. [2] 
While the modes present similar behaviour during cycle 23 up to 
the first half of 2007, the I = and 2 modes behave differently 
compared to the / = 1 modes from June 2007 and onwards. 

The different behaviours observed between values of angular 
degree can be interpreted in terms of the spatial geometry of the 
modes measured with the GOLF instrument. Indeed, Sun-as-a- 
star observations such as the ones made by the GOLF and Mark- 
I instruments are only sensitive to the modes with I + \m\ even. 
This is because the solar rotation axis lies close to the plane of 
the sky from the observing point. The observed / = 1 frequency 
is therefore a weighted measurement of the visible components 
(1=1, \m\ = 1), while in the case of the I = 2 mode, the fit- 
ted frequency corresponds to the weighted measurement of the 
zonal (/ = 2,m = 0) and se ctoral (I = 2,\m\ = 2) components 
dJimenez-Reves et"aU,l2004l) . The sectoral modes are more con- 
centrated along the equator while the zonal modes are most sen- 
sitive to the high latitudes. Thus, the / = 2 modes seen by full- 
disk instruments are more sensitive to the high latitudes of the 
Sun. Moreover, they have larger frequency variations than the 
/ = modes, which average the entire visible solar disk, while 
the full-disk I = 1 mode is more concentrated along the equator. 
The results presented here would then indicate that the process 
responsible for the frequency shifts at high latitudes and related 
to the new solar cycle 24 started late 2007. 
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Fig. 1. Mean frequency shifts (c5vo,i,2) in fiHz of the low-degree solar p modes extracted from the analysis of the 365-day (filled 
circles) and the 91.25-day (plus signs) GOLF spectra. The associated mean error bars are also represented. The corresponding 
10.7-cm radio flux averaged over the same 91.25-day timespan is shown as a proxy of the solar surface activity (solid line). 
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Fig. 2. Mean frequency shifts (6v) in fiHz of the I — 0, 1, and 2 solar p modes (left to right panels) extracted from the analysis of the 
365-day GOLF spectra. The associated mean error bars are also represented. The corresponding 10.7-cm radio flux averaged over 
the same 365-day timespan is shown as a proxy of the solar surface activity (dotted line). 



4. Summary and discussion 

While solar p-mode frequency shifts have been shown to be 
closely correl ated with the solar s urface activity over t h e past 
solar cycles dHowe et all l2002t iJimenez-Reves et aTl [2004; 
IChaminet all 120071) . it appears that during the unpredicted long 
minimum of cycle 23, the frequency shifts of the low-degree p 
modes show unexpected behaviour. Since the second half of the 
year 2007, the frequency shifts are observed to be uncorrected 
with surface activity. While no significant surface activity was 
observed, the frequency shifts of the modes / = and 2 show 
significant increases and those of the I = 1 modes follow the 
extended minimum. The difference between the / = and 2 
modes and the / = 1 modes can be interpreted as a different 
response to the spatial distribution of the solar magnetic field. In 
Sun-as-a-star measurements such as the GOLF observations, the 



modes / = and 2 modes are more sensitive to the high latitudes 
than the / = 1 modes. Such a behaviour of the frequency shifts 
would indicate variations in the magnetic flux at high latitudes 
related to the onset of solar cycle 24. This result qualitatively 
agrees with the m easuremen ts of the cosmogenic isotop es 10 Be 
dBeer et all Il998l) and 14 C dStuiver & Braziunasl Il989l) which 
show an 11 -year periodic variation even during extremely low 
epochs of surface activity, such as the Maunder minimum. 

Also, to estimate in a simple manner when the minimum of 
cycle 23 occurred, we performed a second order polynomial fit 
to the mean frequency shifts (/ = and 2) over the period 2006- 
2009. To do so, we used only independent points from the analy- 
sis of the 91.25-day series. Thus, we obtained that the minimum 
of solar cycle 23 must have happened during the last quarter of 
2007. 
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Table 1. Frequency shifts at each / per unit of change in solar proxy and correlations from the analysis of independent 91.25-day 
GOLF spectra. 
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Fig. 3. Relation of the mean frequency shifts (in j-iHz) between angular degree I for the periods April 1996-June 2007 (open circles) 
and June 2007-May 2009 (filled squares) from the analysis of the 91.25-day GOLF spectra (only independent points were used). 
Left panel: (6vo) against {dvi). Right panel: {5v\) against (6V2). The solid lines correspond to linear fits between the data, and the 
dashed lines represent the 1 : 1 correlation. 



The prediction of the solar cycle properties, such as its 
strength and pattern, give significant l y different conclusions for 
cycle 24 (see e.g., ISvalgaard et al.L 120051: iDikpati & Gilman, 
2006). Most of these models are calibrated using solar observ- 
able, such as the sunspot number, which reflect changes in the 
solar surface. It is clear that there is a time delay between 
the emergence of magnetic fields located probably close to the 
tachocline and the appearance of activity on the surface. The 
physical conditions beneath the surface are thus modified during 
this process. Our findings show that the signature of activity ap- 
pear before in the solar acoustic mode parameters than in other 
solar activity proxies, and therefore must be taken into account 
in the calibration of the solar dynamo models. 

In order to study in more details the spatial dependence 
of the frequency shifts during this peculiar solar minimum, 
this work will be extended to the analysis of the spatially- 
resolved observations from the Global Oscillation Network 
Group (GONG) and Michelson Doppler Imager (MDI) instru- 
ments, which allow the decompo sition into individual m com- 
ponents. iJ^in^z^Re^eTeTal] (|200J) previously showed that the 
values of the frequency shifts for I — 0, 1, and 2 measured in 
MDI data scale to the corresponding spherical harmonic compo- 
nents of the observed line-of-sight surface magnetic field. The 
analysis of the variability of the other p-mode parameters during 
this extended minimum is also currently underway. 
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